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This paper presents a constitutive model and data for materials subjected td large 
strains, high strain rates and high temperatures. The basic form of the model is well suited for 
computations because it uses variables which are readily available in most of the applicable 
computer codes. The materials considered are OFHC copper, Cartridge brass, Nickel 200, Armco 
iron, Carpenter electrical iron, 1006 steel, 2024-T351 aluminum, 7039 aluminum, 4340 steel, S-7 
tool steel, Tungsten alloy and DU*.75Ti (depleted uranium). The data for the material constants 
are obtained from torsion tests over a wide range of strain rates, static tensile tests, dynamic 
Hopkinson bar tensile tests and Hopkinson bar teats at elevated temperatures. The model and 
data are evaluated by comparing computational results with data from cylinder impact tests. 


INTRODUCTION 

During recent years there has been a great deal 
of effort directed at computations for intense impulsive 
loading due to high-velocity impact and explosive deto¬ 
nation. The capabilities of current computer codes have 
been extended to the point that the limiting factor is 
often that of adequately defining material character¬ 
istics for both strength and fracture. A common ap¬ 
proach is to simply repeat the computations with dif¬ 
ferent material characteristics until agreement with ex¬ 
periment is obtained. Clearly it would be desirable to be 
able to characterize materials with a limited number of 
laboratory tests so the initial computations could be 
used with more confidence. This would result in a more 
efficient design cycle and would also give increased in¬ 
sight into understanding the complicated ^ processes 
which ocrur during intense impulsive loading condi¬ 
tions. 

This paper presents a constitutive model which is 
primarily intended for computations. It is recognized 
that more complicated models may indeed giye more ac¬ 
curate descriptions of material behavior. Similarly, vari¬ 
ous models may give better descriptions for various ma¬ 
terials. In many instances, however, the computational 
user cannot readily incorporate complicated and diverse 
models. The result is that a constant ’'dynamic flow 
stress'’ is often used. 

In the following sections test data will be pre¬ 
sented and the approach used to extract the appropriate 
material constants from the data will be explained. This 
will be done in detail for three materials (OFHC copper, 
Armco iron and 4340 steel). Results for the other nine 
materials will be included, but with limited discussion. 
Finally, the model and data will be evaluated by com¬ 
paring computational results to data from cylinder im¬ 
pact tests. 

test data 

The test data are primarily obtained from torsion 
tests over a wide range of strain rates (quasi-static to 


about 400 s" 1 ) and dynamic Hopkinson bar tensile tests 
over a range of temperatures. Static tensile test data are 
also used. The specific torsion testing apparatus is de¬ 
scribed in [1] and a representative Hopkinson bar ap¬ 
paratus is described in [2]. 

Figure 1 shows Hopkinson bar test data at vari¬ 
ous temperatures. The elevated temperatures are ob¬ 
tained by surrounding the in-place test specimen by an 
oven such that the temperatures are applied for several 
minutes prior to testing. Although it is possible to test 
materials to greater strains than those shown in Figure 
1, the Hopkinson bar data cannot be accurately eval¬ 
uated after necking begins to occur in the tensile speci¬ 
mens. Futhermore, at large strains the effects of 
adiabatic heating can also complicate the results. The 
elevated temperatures show a distinct softening effect 
on the strength of the materials. 
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Torsion data for the same three materials are 
shown in Figure 2. Comparable data for the other nine 
materials are presented elsewhere [3,4J. Some desirable 
features of this testing technique are thajt the state of 
stress in the specimen is well defined, large shear 
strains can be achieved without forming geometric in¬ 
stabilities, and a wide range of strain rates can be ob¬ 
tained with the same testing technique. 



AVERAGE SHEAR STRAIN. 




Figure 2. Stress-Strain Data for Torsion Tests at 
Various Strain Rates 


Quasi-static stress-strain data for both tension 
and torsion tests are shown in Figure 3. An equivalent 
tensile flow stress is obtained from the torsion data by 
using the von Mises flow rule; the tensile stress is a = 
and the corresponding tensile strain is e - y/'fS. 




Figure 3. Comparison of Quasi-Static Stress-Strain 
Data for Tension and Torsion Tests 
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The stress for the tension test data is based on the cur¬ 
rent area of the neck, and the strain is defined as In 
(A /A), where A 0 and A represent the initial and current 
areas of the neck. 

At larger strains, after necking has begun, the 
net tensile stress is greater than the tensile flow stress. 
This is due to the presence of hydrostatic tension. The 
equivalent tensile flow stress, obtained from the tension 
data, is approximately determined by using the 
Bridgman correction factor [5j. 

In all three instances the tension data tend to in¬ 
dicate a higher strength than that derived from the tor¬ 
sion data. Therefore, it is desirable to have data from 
both tension and torsion tests such that the discrepancy 
between the two modes of deformation can be identified 
and compensated. 

Figures 4 and 5 show the effect of strain rate for 



Figure 4. Stress vs Strain Rate for the Torsion 
Tests 



Figure 5. Stress vs Strain Rate for the Tension 
Tests 

& 


both torsion and tension. Comparable torsion data for 
the other materials are given in [3] and [4J. In all in¬ 
stances the stress increases as the strain rate increases, 

DATA ANALYSIS AND MODEL DEVELOPMENT 

The model for the von Mises flow stress, a, is ex¬ 
pressed as 

<j « [A + B <s n ] [1 + Cln 6*) [1-T* m ] 

where e is the equivalent plastic strain, c* - i/e 0 is the 
dimensionless plastic strain rate for e 0 = 1.0 s _l and T* 
is the homologous temperature. The five material cons¬ 
tants are A, B, n, C,Am. The expression in the first set of 
brackets gives the stress as a function of strain for e* = 
1.0 and T* ** 0. The expressions in the second and third 
sets of brackets represent the effects of strain rate and 
temperature, respectively. The basic form of the model 
is readily adaptable to most computer codes since it uses 
variables U, «*, T*) which are available in the codes. 


The first step in the process is to determine the 
constants in the first set of brackets. A is the yield 
stress and B and n represent the effects of strain hard¬ 
ening. Since the torsion data include the strain rate of 
interest (e * = 1.0), it is a straightforward procedure to 
obtain the appropriate constants for this strain rate. 
However, these constants are for isothermal conditions 
(T* “ 0) so the effects of adiabatic heating must be con¬ 
sidered in analyzing the torsion data [3,4]. 

The same three constants can also be derived 
from the tension data. Here the starting point is the 
equivalent tensile flow stress from the tension data, as 
shown in Figure 3. This approximate flow stress, ob¬ 
tained by using the Bridgman correction factor, is ad¬ 
justed based on finite element simulations of the tension 
test. The proper flow stress will give computed results 
which agree with the tension test data of Figure 3. Gen¬ 
erally, the Bridgman correction factor gives acceptable 
results, and the adjustments due to the refining com¬ 
putations are not large. 

The flow stress which is eventually determined 
from the quasi-static tension test is for e * - 0.002 and 
must therefore be adjusted for e * - 1.0. This is done by 
increasing the quasi-static values of the constants A and 
B, by the ratio of the stresses at e * = 1.0 and e * = 
0.002, as shown in Figure 5. 

The preceding gives two sets of constants, one for 
torsion and one for tension. For general computational 
purposes, the average values of A, B and n can be used. 
The resulting adiabatic flow stress, based on average 
values of the constants, is generally within about 10 
percent of that obtained by using either the torsion or 
the tension data. 

The effect of thermal softening is shown in Fig¬ 
ure 6. The thermal softening fraction, K^, is simply the 
ratio of the stress at elevated temperature to that at 
room temperature. The data are taken frbm the 
Hopkinson bar tests of Figure 1. The Armco iron is 
'Shown to soften very rapidly with increasing 
temperature, and this is in close agreement with static 
data [3,6]. The OFHC copper shows an approximately 
linear dependence on T*. This is not in good agreement 
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Figure 6. Thermal Softening Fraction vs 

Temperature for the Hopkinson Bar 
Tests 

with static data [3,7], which tend to soften mors rapidly 
as T* increases. There is some evidence, however, that 
the thermal softening fraction tends bo increase as the 
strain rate increases [S]. This trend is consistent with 
the results reported herein. 

The strain rate constant, C, is determined next. 
For the tensile mode of deformation it can be obtained 
from the data of Figure 5. Similarly, for torsion it can be 
obtained from the data of Figure 4. The torsion data also 
allow the combination of the strain rate hardening and 
thermal softening to be evaluated at large strains. This 
is done by numerically simulating the torsion teat at 
high strain rates. Strain hardening, strain rate hard¬ 
ening, thermal softening and heat conduction effects are 


considered [9]. In some instances the strain rate cons- 
tant is adjusted to give better agreement with torsion 
test data at large strains. Figure 7 shows a comparison 
of the computed and test results using the torsional 
strain rate constants. The final strain rate constant is 
taken as an average of those obtained from the tensile 
data and the torsion data. 



Figure 7. Comparison of Test Data and Numerical 
Simulations for the Highest Strain Rate 
Torsion Tests 

The resulting constants for the 12 materials are 
presented in Table 1. The corresponding stress-strain re¬ 
lationships are shown in Figure 8. The temperature for 


Table 1. Constitutive Constants for the Various Materials 


MATERIAL 

GESCXWI0N 

CONSTITUTIVE CONSTANTS FOR 
<r - [a . B« n ][lU^< , ][l - v ' n ] 

fiARONESS 

(ROCKWELL) 

DENSITY 

SPECIFIC 

HEAT 

U/k*KJ 

MELTING 

TEMPERATURE 

(K) 

A 

(MPa) 

B 

(MPa) 

'Tl 

C 

-m 

OFHC COPPER 

F-30 

S9SQ 

383 

1356 

90 

292 

.31 

.025 

1.09 

CARTRIDGE BRASS 

P-67 

' 8620 

its 

1189 

112 

505 

.42 

.009 

1.68 

NICKEL 200 

P-79 

8900 

446 

1726 

163 

648 

.33 

.006 

1,44 

ARMCQ IRON ■ 


7890 

452 

1811 

175 

380 

.52 

.060 

0.55 

CARPENTER ELECTRICAL IRON 

F-S3 

7890 

452 

1811 

290 

339 

.40 

.055 

0.55 

1006 STEEL 

F-94 

7890 

452 

1811 

350 

275 

.36 

.022 

1.00 

2024-T351 ALUMINUM 


— 

875 

775 

265 

426 

.34 

.015 

1,00 

7039 ALUMINUM 

S-7S 

2770 

' 875 

877 

337 

343 

.41 

.010 


4340 STEEL 

C-30 

7830 

477 

1793 

792 

mm 

.26 

.014 

1.03 

S-7 TOOL STEEL 

C-50 

7750 

477 

1763 

1539 

477 

.18 

,012 

1.00 

TUNGSTEN ALL0Y(.07N f. ,03Fc ) 

C-47 

17000 

134 

' 1723 

1506 

177 

.12 

.016 

1.00 

DU-.75Ti 

C-4S 

1*600 

117 

1473 

l 1079 

1120 

.25 

.007 

1.00 
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Figure 8. Isothermal and Adiabatic Stress-Strain Relationships for the Various Materials 













P4.5 


the adiabatic relationships is due to the plastic work of 
deformation. In all cases, the adiabatic stress** reach a 
maximum and then decrease with increasing strain. 
The point of maximum stress is important inasmuch as 
it represents the strain at which localized instabilities 
may begin to occur. 

Some comments should be made about the mate¬ 
rials other than the OFHC copper, Armco iron and 4340 
steel. The Cartridge brass and Nickel 200 have several 
similarities; they are not very sensitive to strain rate. 
[3], the flow stress derived from the tensile data is sig¬ 
nificantly greater than that obtained from the torsion 
data [3], and the thermal softening is leas pronounced 
than for the other materials. The thermal softening data 
obtained from the Hopkinson bar tests gives m « 1,68 
(Cartridge brass) and m 35 1.44 (Nickel 2Q0X This is con* 
sistent with static data [3,10, 11]. 

No Hopkinson bar tests were performed on the 
Carpenter electrical iron and it is assumed that its 
thermal softening characteristics are similar to that of 
the Armco iron. Hopkinson bar tests were performed on 
the 1006 steel and some very unexpected results were 
obtained. At lower temperatures (0 < T* < .2) tile 
thermal softening data are similar to that of the Armco 
iron. At higher temperatures, however, (.2 < T* < .4) 
K T increases sharply as the temperature increase*. Sin¬ 
ce these results cannot be adequately explained, the 
thermal softening fraction is assumed to vary in a linear 
manner (m = 1.00). The other constants for the 1008 
steel are obtained from the torsion dkta only [3]. 

The data for the remaining material* (2024-T351 
aluminum, 7039 aluminum, S-7 tool ateel, Tung*ten al¬ 
loy and DU-.75Ti) are also obtained from torsion data 
only [4]. No tensile data are included. The thermal sof¬ 
tening fraction is again assumed to vaiy in a linear 
manner. 

EVALUATION WITH CYLINDER IMPACT TESTS 

An evaluation of the model and data can be made 
by comparing computed results with data from cylinder 
impact tests. These t^sts provide strain rate* in excesa 
of 10 s s" 1 and strains in excess of 2.0. Also, a range of 
-tensile and shear deformations is experienced. 

The results of these comparison* are shown in 
Figure 9. The computations were performed with the 
EPIC-2 code. For these computations an average pre*- 
sure was used for each set of adjacent triangular ele¬ 
ments such that the excessive stiffness sometime* as¬ 
sociated with triangular elements is essentially 
eliminated [12]. The Armco iron and 4340 *t*el show 
very good agreement. The OFHC copper results do not 
agree as well, but even these are acceptable for moat ap¬ 
plications. It should be emphasized that th* result* of 
cylinder impact tests have not been incorporated into 
the data, and they represent a totally independent check 
case. Photographs of the impacted cylinder* are shown 
in Figure 10. 

SUMMARY and conclusions 

A computational model for flow stress has been 
presented which includes the effects of strain hardening, 
strain rate hardening and thermal softening. Some well- 
defined laboratory tests have been performed on 12 ma- 




NOTEJJ: 

-INITIAL LENGTHS L„«2.SACM UNLESS OTHERWISE NOTED 

-INITIAL DIAMETERS D*» 7B2CM 

-ELASTIC STRAIN CONTOURS SHOWN AT INTERVALS Of 0.5 

- TEST DATA OENOTED BY DOTS { • • • ) 

Figure 9. Comparison of Computed Shapes and 

Test Results for Cylinder Impact Tests at 
Various Velocities 

terials to obtain da‘ta for this model. An independent 
evaluation has shown the model and data to give good 
computational results for a range of cylinder impact con¬ 
ditions. 
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